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Abstract 

We study the transition amplitudes in the state-sum models of quantum gravity in 
D = 2,3,4 spacetime dimensions by using the field theory over formulation, 
where G is the relevant Lie group. By promoting the group theory Fourier modes into 
creation and annihilation operators we construct a Fock space for the quantum field 
theory whose Feynman diagrams give the transition amplitudes. By making products 
of the Fourier modes we construct operators and states representing the spin networks 
associated to triangulations of spatial boundaries of a triangulated spacetime mani- 
fold. The corresponding spin network amplitudes give the state-sum amplitudes for 
triangulated manifolds with boundaries. We also show that one can introduce a dis- 
crete time evolution operator, where the time is given by the number of D-simplices in 
a triangulation, or equivalently by the number of the vertices of the Feynman diagram. 
The corresponding transition amplitude is a finite sum of Feynman diagrams, and in 
this way one avoids the problem of infinite amplitudes caused by summing over all 
possible triangulations. 
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1 Introduction 



The idea of constructing a quantum theory of gravity by assuming that the spacetime is 
discrete at the Planck scale is appealing because it resolves by definition the problem 
of UV infinities of the conventional quantum field theory approach. However, the 
development of this idea was hindered by the complexity of the obvious candidate for 
such a theory, i.e. Regge simplicial formulation of general relativity (GR) (jH]. 

In the past decade a new class of simplical gravity models have been developed, 
which are based on the representation theory of the relevant symmetry group G. This 
group is SO{D) for the Euclidean gravity case or SO{D — 1, 1) for the Lorentzian 
case. The prototype model was the Ponzano-Regge model of 3d gravity [|^], which 
was made mathematically well defined in the work of Turraev and Viro by passing 
from the 50(3) to the quantum SU{2) group [|]]. This was generalized to 4d in the 
approach of topological state sum models [Q], which is a mathematical formalism for 
constructing topological invariants for manifolds out of colored triangulations, based 
on the category theory. In physics terms, this is a way of constructing partition func- 
tions for topological field theories. The state-sum formalism was based on the papers 
of Boulatov and Ooguri, who realized that the topologically invariant state sums can be 
generated as Feynman diagrams of a field theory over [|5|, ^. Their approach was 
motivated by the result from string theory that the partition function of a discretized 
string theory can be represented as a sum of Feynman diagrams of a matrix field theory 
(for a review and references see [^). 

Based on these developments. Barret and Crane have proposed an approach for 
obtaining the partition function for general relativity (GR) in D = 4, which is a non- 
topological theory, from a state sum of a topological theory, the BF theory []|, ^. 
The BC model represents a very interesting formulation of discrete quantum gravity, 
since the local symmetry group of the spacetime plays a crucial role. One colors 
a triangulation T of the spacetime manifold M with a certain class of irreducible 
representations (irreps) of the Lorentz group, and the partition function for a given 
T is the sum over the amplitudes for all colorings. Given the partition function, one 
can compute the transition amplitude associated with a triangulation of M with spatial 
boundaries which are colored with a fixed set of irreps, by summing over the internal 
irreps, exactly as in the topological cases [||, |m pT| ]. 

Since a field theory formulation of the BC model exists [|T2|], such an amplitude 
would correspond to a Feynman diagram with external legs. In the quantum field the- 
ory (QFT) formalism such diagrams contribute to the matrix elements of an evolution 
operator, the S-matrix. This evolution operator acts in the Hilbert space of the QFT, 
which is the Fock space constructed out of the creation and the annihilation operators. 
The aim of this paper is to develop these concepts for the BC model. 

A further motivation is to construct the states which correspond to the colored tri- 
angulations of the spatial boundaries, i.e. the spin network states. This is based on 
the property of the QFT formalism that one can associate states to operators via the 
Fock space vacuum state, so that one could in principle construct the operators corre- 
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spending to spin networks. This would be also the first step for the third quantization 
formulation, i.e. a quantum theory where the number of the universes is not fixed. 

Furthermore, one can define the transition amplitudes between these spin net states 
by summing over the appropriate Feynman diagrams with external legs. Since such a 
Feynman diagram (FD) corresponds to a triangulation of the spacetime with bound- 
aries, there are several amplitudes one can associate to a transition from a state on the 
initial surface to a state on the final surface. The standard amplitude is given by the 
sum over all possible triangulations, or equivalently over all FD, but we argue here that 
it makes sense to consider the perturbative part of that amplitude, which is given by 
the sum of the FD with n vertices, or triangulations with n simplices, as the transition 
amplitude for the time interval of n Planck units of time. In this way one introduces a 
discrete time variable, and also regularizes the amplitude, since one avoids the infinite 
sum over all triangulations, which is a generically divergent expression. 

In sections 2,3 and 4 we study the simplical field theory for a 2d discretized gravity 
theory based on the 2d BF theory, since it is a good toy model for illustrating and 
developing our ideas. In section 5 we further develop and explore the general concepts 
introduced in the case of the 2d model, where we consider the 3d simplical field theory 
of Boulatov. In sections 6 and 7 we apply these concepts to the field theory formulation 
of the EC model. In section 8 we present our conclusions. 

2 D=2 model 

We start with the D = 2 model, since it is very useful for understanding and illustrating 
the basic concepts. The D-dimensional BF theory is given by the action 

Sbf = [ <BAF> (1) 

where B is a (D — 2)-form,i.e. scalar field in 2d, and F = dA + AAAis the curvature 
two-form for the connection one-form A. A and B take values in the Lie algebra of 
the Lie group G and < , > is the corresponding invariant bilinear form. 

For the case of 2d gravity the relevant groups are the 2d Poincare group IS0{1, 1), 
then the 2d anti-de-Sitter group 5*0(1, 2) and its Euclidean versions 5*0(3) and SU{2). 
The corresponding BF theory gives the Jackiw-Taitelboim theory which follows 
from the identifications 

A = ujJo + e^J± , B = B°Jo + B^J± (2) 

where to is the spin connection, are the zweibeins, B^ is the dilaton, and J are the 
Lie algebra generators, so that 

Sbf= [ B^dco + Xce^ Ae-) + B^{de^±ue^), (3) 
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where Ac is the cosmological constant. Hence the enforce the zero-torsion condi- 
tions which give uj = uj(e) and one ends up with the Jackiw-Teitelboim action. We 
will consider the compact G case, i.e. Euclidean 2d metrics, since the formulas are 
simpler, and the basic ideas are the same. 

The field theory which generates the partition functions for the triangularizations 
of the two-manifold M as Feynman diagrams is given by 

S2 = \ I d'^9^'^i9i,92) + ^ [ d^gcpigi, 92)^(92, 93)^(93,91) , (4) 
^ Jg^ -^i Jq3 

where satisfies 

^(91,92) = V{9i9,929) , (5) 
and A is the perturbation theory expansion parameter. Hence 

f{9i, 92) = f{9i ■ 92^) = Yl ■ 92^) , (6) 

A,a,l3 

where the last formula follows from the Peter- Weyl theorem. A denotes an unitary 
irreducible representation of G, 1 < a,P < dj^ = dim A. The Fourier coefficients 
(p^ will be important for the construction of the spin-network states and transition 
amplitudes, and they are the analogs of the creation and the annihilation operators 
from the particle field theory. 

By inserting (H) into (@) and by using the orthonormality relations 

/ d9{DU9))*Da'A9) = ^S^'^'5a,a'S,,,^ , (7) 
Jg '^a 

and the complex conjugation relations^] 

iDU9))* = i-^)''-''-'D^a-pi9), (8) 
we obtain for the kinetic part of S 



A,a,l3 A,a,f3 



(9) 



while for the interaction part we get 

3! 



= ^ E df<P%<Pl<P^^ . (10) 



A,o,/3,7 

Note that (|]) induces the following reality condition for the Fourier components 

= , (11) 



tXhis form is for the SU{2) case, so that A = 2j = 0, 1, 2, ... and c^a = A + 1. 
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which is a non-AbeUan generaUzation of the reaUty condition for the usual Fourier 



modes: a* = a^p. 



As far as the Feynmanology is concerned, it will be useful to consider a general 0^ 
theory given by the action 

-S" = ^ ^ (t>i(t>jCij + ^ ^ (t>i(t>j(t>KVijK , (12) 

I,J ' I,J,K 

where the indices /, J, K belong to a general label set. 

The perturbative expansion an the Feynman diagrams are generated by using the 
generating functional 

exp(i(5[0]+ J] , (13) 

I 

where Vcf) — Y\_i d(f)i. One can write 

Z{J) = J2Ji,---Jj^Gj,...j^ (14) 

where the Green functions G can be written as 

n\ ^-^ ml J 

m 

ES«s.r0,r--^o = SEG&. ■ (15) 



m — ' m\ ni 

m m 



The perturbative components G^"^^ can be calculated by a repeated differentiation 
of the Gaussian integral 

Zo[J] = J ^^0exp = Aexp (^-i^J^JjCyjjj^ (16) 

with respect to the currents J, where A is a constant proportional to det " 2 1 1 C*^ j 1 1 . One 
can then show that the result for G^'^^ can be expressed by using the classical analog 
of Wick's theorem from quantum field theory, i.e. one makes all possible pairings in 
the string of cp's given by 6*1, ■ ■ ■ SyCpi-^ ■ ■ ■ 0/„ and then assigns to each pairing (0/, (pj) 
a number equal to Cj'j, and then takes the product of these numbers and sums over all 
possible sets of pairings. 

Wick's theorem is the basis for the Feynman diagrammatic expansion since topo- 
logically non-equivalent sets of pairings can be denoted as trivalent graphs, which are 
the Feynman diagrams. For example, the term in G/j^^ given by 

^ijKL — # ^ GjpGj],Gj^\,Gi}^,Gj2n^i'J'mVk'L'n (17) 

r,J',K',L',M,N 
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where # is the number of equivalent pairings, will correspond to the H diagram. In 
order to obtain its value, which is the corresponding S -matrix contribution, or the 
transition amplitude, one has to remove the external legs propagator factors. This 
corresponds to calculating 



AljKL — CiiiCjj'CkK'ClL'^I'J'K'L' 
I'J'K'L' 

= # Cj^f^VijuVKLN ■ (18) 
MN 

The partition function is given by (e*'^" ) and it is equal to the sum of all vacuum 
diagrams. The connected diagrams are generated by W[J'] = logZ[J']. If W = 
J2n W^n-^"' then the connected vacuum diagrams are given by 

W2 = (theta) + {dumbbell) 

W4 = {Mercedes — Benz) + {cylinder) + ■ • • 

Wq = {prism) + ■ ■ ■ , (19) 

where 

{theta) = WijkVlmnC'^C'^'C''^ 
{dumbbell) = WijkVlmnC^^C^^C^'^ 
{Merc. -Benz) = WijkVlmnVqpqVrstC'''' C''''C'^''C'''C''C^\{2Qi) 

and so on, where C^"^ = Cjj and 7^ is the corresponding combinatorial factor. 

Hence the vacuum Feynman diagrams denote the rules of forming invariants by 
contracting the indices in the product of tensor quantities V, i.e. given a product 
^/i JiA'i • ■ • VinJnKn^ the corresponding Feynman diagrams represent all topologically 
non-equivalent ways of contracting the indices in this product with C^'^ = Cjj. The 
matrix C^'-' is the propagator in the field theory language. If a set of indices {/, J, ...} 
remains uncontracted, then this gives a Feynman diagram with external legs, which 
contributes to the transition amplitude for a set of states carrying the quantum num- 
bers {/, J, ...}. 

Consider the {n + m) -point Green function 



jn+m 

(n + m)! 



= , , .M '^h, ■ ■ ■ hjj, ■ ■ ■ <PjJ ■ (21) 



It will generate all FD with (n + m) external legs. Note that the subclass of these 
FD containing all connected FD which have a string]] of n legs and a string of m legs 
will contribute to the amplitude for a transition from a state labeled with quantum 
numbers Ji, /„ to the state labeled with quantum numbers Ji, Jm- This can be 



■^By a string we mean that there is no vertex with an internal Hne between the legs of a string. 
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explicitly realized by using the concept of creation and annihilation operators from 
the operator quantization formalism. The index label set / can be always split as 
I — j, —j, k,l} such that 

0- = 0_i , 0* = , (t>l^(t>k , 0* = -(t>i ■ (22) 

Then promote 0±j and 0±j modes into operators 0±i and \/— T0±j satisfying 

ct>\^U , (23) 

where i.j G {i, j}. The zero-modes 0^ and 0; promote into operators |(0ifc + ^l) and 
1(0/ - 0|) where 

[</'^,</'r]=<^M ' ' Wl&iKl}, (24) 
and all other commutators are zero. 

These operators naturally act on the Hilbert space of states 

4^---4|0) , 0HO) = O, (25) 

where now i' e {i, j, A;, /}. Then the amplitude for a transition from the initial state 

= |0) (26) 

to the final state 

\^2) = <py--<p], |0) , (27) 

is given by 

Ai2 = (0| ■ ■ ■ <P,>^ cxp(z : :)0|, • • ■ ^ |0) , (28) 

where : 5*^, : is the normal-ordered operator Sy with respect to the vacuum |0), and all 
the subsequent manipulations are the same as in the particle field theory case, i.e. one 
uses the operator formulation of Wick's theorem. 

Note that one does not need to go to the operator formalism in order to compute the 
amplitudes. One can simply say that there are abstract states | Ji...J„) and the transition 
amplitude from the state | Ji ■ ■ ■ f^i) to the state IL2 ■ ■ ■ M2) will be given by the sum 
of the corresponding Feynman diagrams. We will denote the result as 

Au = (L2 • • • M2 1 eMiSv) \Ji---Ki) . (29) 

We will also use the notation 

^12 = (0L2 • • • 0M2 exp{iSy)(f)j^ ■■■(pKi) , (30) 

because it indicates that the transition amplitude is obtained from the Green's function 
by taking only the contractions corresponding to connected Feynman diagrams with a 
string of n and a string of m external legs, which we will denote as (n, m) FD. 
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The expression ( j^ ) is calculated perturbatively in A, where at order n one has to 
calculate 

4? = (0L.---<^M.(^.)>Ji---</)i^i) . (31) 

In the case of field theories of particles, the physically relevant quantity is (P9|), which 
gives the transition amplitude from a state in a distant (infinite) past to a state in a 
distant (infinite) future. Therefore in order to obtain meaningful amplitudes one has 
to sum all Feynman diagrams up to a given order of perturbation theory. However, in 
the case of field theories associated to simplical gravity, we will argue that even the 
individual Feynman diagrams, or certain sub-classes of diagrams, can have a physical 
interpretation. This means that it is not necessary to sum over all possible triangula- 
tions interpolating between the initial and the final surface. 

Note for example that in the case of the vacuum diagrams, a single diagram V is 
the partition function for the simplical complex T where F = T^* is the dual one- 
complex to T. If one thinks about the simplical gravity as the fundamental theory,i.e. 
not as an auxiliary device for defining the path integral of the continuum theory, then 
each equivalence class of simplical complexes could have a physical meaning. This 
is trivially satisfied in the case of topological theories, where due to the triangulation 
independence, the continuum theory is the same as the simplical theory, and hence the 
partition function for a single simplical complex T(M) is the partition function for M. 



3 Partition function 

In order to develop these ideas further, we examine in some detail our 2d simplical 
gravity model. Consider the Mercedes-Benz (MB) diagram. It is a dual one-complex 
for a tetrahedron t, if the t is considered as a 2d simplical complex. In that case 
t = T{S'^), and hence the value of the MB diagram will be the partition function for 
a sphere 5*^. The kinetic term @ implies that the propagator will be given by two 
parallel lines, representing the basic contraction 

{<P% = (-l)^-"-''rfA5^'^'5„,-„'5A-/5' . (32) 

The interaction term (|T^ implies that the vertex will be given by a triangle whose 
corners are cut, where the indices G Aj are assigned to each edge. One then joins 
the vertices by the propagator lines, and in this way one obtains a ribbon MB graph. 
The value of the diagram will be given by the product of the edge delta functions 
multiplied by the propagator factors and the vertex factors. In general there will be 
2E delta functions, which will be multiplied by a factor (—l)^^{d\)^~^^ , where E is 
the number of edges and V is the number of vertices of F. For the MB case E = 6 
and = 4. By summing over the indices along a loop of F, the corresponding delta 
functions give the factor c/a, so that 

Z{T) = ^(dA)^-2^+^(-l)^^ 

A 
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A A 

where F is the number of loops in the ribbon graph T (F{MB) = 4), which is the 
same as the number of faces in T*, and we have used that 2E = 3V for trivalent 
graphs. 

Since x(r) = = x(^) = 2 — 2g, the Euler number, is a topological 

invariant, where g is the genus of M, Z would be a topological invariant provided that 
there is no sign factor (-1)^^. It comes from the sign factor in the propagator (521). 
There are several ways one can get rid off this sign factor. One can change the kinetic 
term into 

= ^ E , (34) 

A,a,l3 

SO that the sign factor is removed from the propagator. Although this choice gives 
the correct answer, this action does not follow in general from the integration over 
the group, unless one has a group where one can choose a basis where {D^p)*{g) = 
D^^_p{g). Similarly, one can choose 

^'^ = ^E^a' (35) 

and this would follow from the group integration if the group G allows a basis where 
the matrices D{g) are real. In both cases we will write the propagator as 

(0(12) 0(34)) = 5(13)5(24) . (36) 

One can avoid the propagator factor ci]^ by rescaling 0, which then changes the 
vertex factor io d^^'^ . Either way the propagator ( |36| ) gives the required result 

Z{M) = ^{dj.Y^'''^ . (37) 

A 

If one wants to have an integral over the group expression which is valid for arbi- 
trary G, then one can choose 



^k = \j <^{9i, 92)^(92, gi) = ' 

A,a,/3 



tp^i. ■ (38) 



This choice gives a twisted propagator 

(0(12) 0(34)) = 5(14)5(23) . (39) 

The difference between the twisted propagator case ( P9| ) and the untwisted case ( P5[ ) 
is that the same graph will not describe the same surface in each case. For example, in 
the case of the theta graph, the untwisted propagator gives x = 2, a sphere, while the 
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twisted propagator gives x = 0, a torus. In the untwisted case the torus is obtained by 
braiding two edges of the theta graph, which for the twisted case gives a sphere. 

Note that one can impose the symmetry condition Lp(l, 2) = (^(2, 1), which to- 
gether with the reality condition implies 

(0a/3)* = </>i , (40) 

i.e. (p^ are hermitian matrices. This gives a group theory generalization of the hermi- 
tian matrix models used in string theory. 



4 Transition amplitudes and spin networks 

Now let us try to formulate a transition amplitude from a spatial section Ei to a spatial 
section S2, which are boundaries of M, i.e. we cut two holes in M. Consider the 
case when the corresponding simplical complex is a prism. Then the base triangles 
correspond to triangulations of Sj = S^, where is a circle. Each side of the prism is 
divided by a diagonal into two triangles, and this defines our triangulation T{M) where 
M = S"^ a two-sphere. The corresponding Feynman diagram is given by the dual 
graph, which is a circle with six legs. Its value can be obtained by making appropriate 
contractions in 

• (41) 

If we take ( pT] ) as it is, we will obtain an expression which is not G invariant, 
because of the uncontracted representation indices. In order to remedy this, consider 
the following quantity, or the operator 

$3=/ d^gip{gi,g2)(p{g2, 93)^(93,91) ■ (42) 
It is proportional to Sy, so if we use ([T0|), we get 

$3=5^ t/XVV^,0t = E^A''^3(A) . (43) 

A,Q,/3,7 A 

We then propose that $3(A) be an operator associated to a Id spin network given by a 
triangle whose each side is colored by the irrep A. Note that $3 (A) is defined up to a 
scale factor, so we can write 

*^'''>=]V^ 5: = ]^rr(^-)» . (44) 

A natural normalization factor would be A^(A) = cf/'^. 
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One can then associate states |$3(A)) to operators $3 (A), and the corresponding 
amplitude will be given by 

4? = (<^'3(A)| |$3(A)) = (<f3(A)(5.)e$3(A)) 

= ]v4^E('^^l/3,0^,,,0^...(^.)V^:,,0^,7>^,.,) . (45) 

If we think of \E'3 as a quantity associated to a triangle whose each vertex carries a 
group element g, then $3 (A) can be interpreted as a quantity corresponding to a spin 
net formed by a spatially dual diagram to the triangle, where the duality is with respect 
to a Id space, and hence one obtains a triangle as a dual complex. Note that we can 
think of ip{gi,g2) = v^(l,2) as a quantity associated to the edge at whose ends are 
attached the group elements gi and g2. Therefore the quantity 

/ d-g^{l,2M2,3)---^{n,l) (46) 

will correspond to an n-sided polygon, whose each vertex carries a group element. 
Hence $„ will generate a spin network variable $n(A) oc Tr (0^)", which will corre- 
spond to a closed Id spin network with n edges colored by A. 

It is easy to see that the non-zero transition amplitudes will come from the pertur- 
bative amplitudes 

= {^mW\ (S.r = ($™(A)(5.)'=<f/(A)) . (47) 



5 Time variables 

Given this setup, and our remarks at the end of section 2, we can say something about 
possible time intervals for the amplitudes. The amplitude for the prism ( p5| ) could 
be interpreted as a transition amplitude from the state $1 to the state $2 in one unit 
of discrete time, so that the transition amplitude in n units of time would be given 
by a sum of FD corresponding to triangulations consisting of n prisms put on top of 
each other. This suggests a discrete time variable T to be the minimal distance, i.e. the 
number of edges, from the initial to the final triangle in the corresponding triangulation. 
However, this time variable will make sense only for triangulations where the minimal 
distance between any two points of the initial and the final triangle does not exceed 
n + 1 . In general case these are the triangulations where for every point p of the inital 
polygon P we have 

min {d{p,p') \p' E P'} = n , (48) 

where P' is the final polygon. 

Note that the triangulations which satisfy (^Sj) can be sliced into (n — 1) polygons 



Pfc, such that the distance ( p8[ ) is one for each pair (P^, Pk+i), where Pq = P and 



Pn = P' . It is easy to see that there will be infinitely many such triangulations for 
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every n > 2. Hence if one tries to define the amplitude A{n) for n units of time as 
a sum over triangulations with n — \ slices, then for n > 2 the amplitude A{n) will 
contain infinitely many FD. Therefore one will have the same problem as in the case 
of the unrestricted amplitude 

Az^ = ($„(A)exp(25,)$,(A)) . (49) 

One can then try to restrict the infinite sum of FD in A{n) by using the perturbation 
theory, i.e. by taking only the FD up to a given order of A. However, one has to insure 
that the corresponding amplitudes are consistent. Namely, the amplitudes have to obey 
the composition rules, which come from the procedure of pasting two cylinders. This 
is equivalent to requiring that the amplitudes are given by the matrix elements of an 
evolution operator U (T), which satisfies 

f/(Ti)[/(T2) = f/(Ti+T2) , (50) 

where T is the evolution parameter, i.e. the time variable. 

For example, in the case of particle field theories, one has U{T) 
exp(iTif^), where the split Sy = THy is possible because of the fixed spacetime 
background structure. In the case of simplical field theories, this split is not natural, 
and it is not clear how to do it. One can then use the minimal distance time variable 
T, and the corresponding U (T) would be defined by requiring that its matrix elements 
are given by the restricted amplitudes A{n). However, it is not clear why would such 
an operator satisfy (^D|), and it is difficult to check this. 

On the other hand, note that U(T) = (S^)'^ satisfies the composition law (pO|). 
Hence the amplitude ( P7[ ) can be interpreted as the transition amplitude for T = k 
units of time, where the time variable is given by the number of vertices in the corre- 
sponding Feynman diagrams, or equivalently, by the number of triangles in the corre- 
sponding triangulations. Since the amplitude ( P7[ ) is given by a finite sum of FD, the 
new evolution operator will give finite amplitudes for finite time intervals, provided 
that the individual FD are finite. This can be achieved by going to the quantum group 
formulation. 

Since this choice of time is related to the covariant quantization procedure we are 
using, i.e. we are not splitting the spacetime into space and time, as in the canonical 
quantization, we will call this evolution parameter the covariant time. 

However, there is a peculiarity with the covariant evolution operator. Since it is 
natural to define Sy to be a hermitian operator in the field theory Fock space, so that 
e*'^" is a unitary operator, then U (n) = (5*1,)" is going to be only a hermitian operator, 
but not a unitary operator. This would mean a non-unitary time evolution of states. 
However, one has to keep in mind that we are not dealing here with the usual domain of 
quantum mechanics, i.e. objects in a fixed background spacetime. We have here states 
describing the whole universe, which is in our toy model represented by a colored 
triangulation of S^, and hence one cannot apriori expect that the same rules apply as 
in the standard quantum mechanics. 
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6 D=3 model 



We now examine our constructions for the case of 3d gravity. The partition function 
for the 3d BF theory is generated by the FD of a (j)"^ field theory on whose 
action is given by 

Ss = ^ [ d'g v\l,2,3) + ^ [ rf^^7<^(l,2,3V(l,4,5)<^(2,5,6V(3,6,4) , 

(51) 

where we use a short-hand (p{l, 2, 3) = ^{gi, g2, g-s.) and satisfies 

^{9i,92,g3) = ^{gig,g2g,g3g) ■ (52) 

The condition (Q) implies 

^{91,92,93)= / dg<p{gig,g2g,g^g) , (53) 



JG 

where is an unconstrained field. 

Note that the form of S3 allows the following (category theory) interpretation. If 
we associate ^9(1, 2, 3) to a triangle whose edges are labeled with the group elements 
gi, g2, g3, the interaction term can be associated to a tetrahedron whose edges are 
labeled by the group elements gi, ...,ge, and the triangles have orientations provided by 
the outer normals. With this interpretation it is natural to require the cyclic symmetry 
for (p, i.e. 

<^(1,2,3) = V9(2,3,1) = <^(3,1,2) . (54) 
However, as in the 2d case, one does not have to impose this symmetry. 



The Fourier modes expansion of ( p3| ) gives 



V = Y. V^A.^A.^Aa <Ptlt!a!Dt±a!{9l,92, 93) , (55) 
A, a 

where 

I A1A2A3 _ j,AiA2A3 ^AiA2A3 fcf-^ 

raia2a3 — V^ai/3ia2/32a3/33^/3i/32/33 ' ^-^^^ 

Dtla!t!i9u 92, 93) = E <V(^iX% (^2X^,3 (^73)C^%tft. (57) 

The coefficients C are basis components of the intertwiner map 63 

63 : V{Ai) ® V(A2) ® l^(A3) V{0) , (58) 
where V{0) is the subspace of the singlets, so that 

|0; Ai, A2, A3) = J2 ® W2) ® las) , (59) 
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and \a) are vectors of an orthonormal basis. 



In order to obtain ( p5| ) one needs the relation 

dgU''^ig)(E)U'"^ig)(E)U'''^ig) = tsi3^ , (60) 



G 



which in the basis (|59|) becomes 



/ d9Dtl,,{9)Dtl,S9)Dtl,,{9)=C^^^^^^^ • (61) 

Jg 



By inserting ( p5| ) into 5*3 we obtain 

2 



~ 2 ^ |^a(123) 



,A(123) 



A, a 

and 



C _ A Vr-n2:.(j.-".)^A(123) ,A(145) ,A(256) ,A(364) JAi A2 A3 

- 41 ^) V>a(123) Va(-U5) ^a(-2-56) <Pa(-3-6-4) ] ^5 Ag 

A.a ^ 



where we have used a shorthand X{jkl) = XjXkXi and 



Al A2 A3 
A4 A5 Ae 



is the 6j symbol. We also denote it as 6j(Ai ■ ■ • Ag). 
Since the reality condition can be written as 



then the kinetic term ( p2\ ) gives for the propagator 



(62) 



(63) 



\^r-nS,0-ft)r^(123) ^A(145) ^A(256) ^A(364) 
/ A ^) ^13(123) '^[^(-145) ^l3{-2-56) ^ l3{-3-6-4.) v"'*; 

/3 



«;s)'=(-i)=*--'C?v„. (65) 



Note that the 3d propagator contains a sign factor analogous to the one in the 2d case. 
Since for our purposes this sign factor, as well as the related sign factors, will not be 
important, we will omit them from the formulas. We will then write the formulas as 
for the real D{g) case. 

The Feynman diagrams will be given by the four-valent graphs F, where each edge 
is represented by 3 parallel lines, while each vertex is represented by a tetrahedron 
whose corners are cut and the corresponding edges are joined with the propagator 



lines. The Feynman rules which follow from ( |6^ ) and ( |63| ) consist of assigning an 



index ctj G Aj to each edge of each tetrahedron, as well as the corresponding 6j 
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symbol. One then assigns a delta function for each propagator line and multiplies 
all the delta functions with the 6j symbols. If T has \V\ vertices and \E\ edges, one 
then obtains a product of 6j symbols and 3\E\ delta functions. By summing over 
the indices for each loop| F of F one obtains a d^Ap) factor, so that the value of the 
diagram is given by 

^(r) = E n ^(^^) n 6J i^F.iv) ■ ■ ■ AF.(y)) . (67) 
Ap Per Ver 

For example, consider a vacuum diagram given by the pentagram graph. The pen- 
tagram is a dual one-complex to the 4-simplex if the 4-simplex is considered as a 3d 
simplical complex. This 4-simplex can be considered as a triangulation of the space- 
time manifold S^. 

If the graph F is dual to a triangulation T(M), we can rewrite Z as 

ziM) = E n ^(^^) n 6j (Kit) ■ ■ ■ Kit)) (68) 

Ae eer ter 

where e are the edges of T, which correspond to the faces F of F. t are the tetrahe- 
drons of T, which correspond to the vertices V of F. In this way one reproduces the 
Ponzano-Regge formula, where the edges e of the triangulation T are colored with the 
SU (2) irreps, and the 6j symbols are assigned to the tetrahedrons t. The sum over the 
colorings is divergent, and it can be regularized by using the quantum group SUq{2) 
[|^]. The result is triangulation independent, so that Zq is a topological invariant of M. 

As far as the transition amplitudes are concerned, they will be given by a straight- 
forward generalization of the D = 2 amplitudes. For example, if we wish to calculate 
the transition amplitude from an initial surface Si = 5^ to a final surface S2 = S^, 
we take a 3d simplical complex with two disjoint boundary tetrahedrons ti, which are 
taken to be the triangulations of the initial and the final surface. The corresponding 
Feynman diagram will be given by the dual one-complex, i.e. a four-valent graph with 
8 external legs. For the invariant amplitude one has to consider the D = 3 analog of 



(|42|), i.e. to a tetrahedron whose edges are labeled with group elements we assign 



$4= / A<^(1,2,3)<^(1,4,5)<^(2,5,6)<^(3,6,4) . (69) 
By expanding y^'s in the Fourier modes we obtain 

^ _ , A(123) , A(145) , A(256) , A(364) f Al A2 A3 1 

^4 - 2^ 0«(123) V^a(145) V>a(256) V>a(364) i A, A. A. ( ' 



^a(123) V^a(145) V^a(256) V^a(364) ^ 
A, a ^ 



Hence we define 



^ r A A ^ - ^ \^ J,A(123) MU5) , A(256) ,A(364) 

^a[Ai, Aq) - <Pa(123) <Pa{U5) <Pa{256) V^a{364) ' 



(71) 



§This is also a face if F is considered as a two-complex. This two-complex is called a spin foam 
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where a natural normalization would be A^(A) = ^yd{Ai) ■ ■ ■ (^(Ae). 

We then associate the quantity $(Ai, ■ ■ ■ , Ag) to the spin net corresponding to a tri- 
angulation of 5^ given by a tetrahedron t. The spin net graph will be given by the dual 
one-complex to t, where t is considered as a 2d simplical complex. One then obtains 
the Mercedes-Benz graph, and its edges will be labeled with the irreps Aj. Hence a 
perturbative amplitude for the transition from |$(Ai, ■ ■ ■ , Ag)) to |$(A'^, ■ ■ ■ , Ag)) will 
be given by 

Ajin) = ($(A;, • • ■ A[j)(5,)"<l>(Ai, ■ ■ ■ , Ag)) , (72) 

where n > 8. 

If we consider a T(S'^) which is two tetrahedrons with a common face, then the 
dual graph is a prism. If we attach to the edges of T the group elements gi, ■ ■ ■ , gg, we 
can then construct an invariant operator 

^g=/" A</'(l,4,2)^(2,5,3)¥p(l,3,6)(/.(4,7,8)</.(5,8,9)(^(7,6,9) . (73) 

By expanding into the Fourier modes we obtain $g(Ai ■ • • Ag), which is an operator 
associated to the spin net given by a prism whose edges are labeled by the irreps 
Ai • ■ ■ Ag. In this way we obtain operators and states for 2d spin nets whose graphs are 
dual to triangulations of a boundary surface S. This surface represents a space part of 
the spacetime manifold M. 

Therefore we denote the $„(A) operators and states as $(7, A), where 7 denotes 
the spin net graph, which is dual to T(S), and A denotes the labeling of the edges of 
7 with the unitary irreps of G. One can write 

^(ry A) - ^ '^TT J)A(eiMe2We3(i»)) .ja. 
^ ~ iV(A) ^-l-l*'^«(eiWe2He3W) ' ^'^^ 

where v are the vertices of 7, and ei{v) are the three edges coming out of the vertex v. 

The amplitude ( [72| ) is equal to a sum of FD with n vertices and 6 + 6 external legs. 
The value of the each diagram in that sum can be written as 



AiT) = NiT) J2 6j(A; ■ ■ ■ A'g) J] ^(A.) J] 6j(Ae-,(t) ■ ■ ■ A,,(,-))6j(Ai ■ ■ ■ Ae 
where 



(75) 



= iV(A;...A-)Ar(A,...Ag) = ^ d{A,) ■ ■ ■ d{A,)d{A,) ■ ■ ■ d{A,) , (76) 

and t are the interior tetrahedra and e are the interior edges of a triangulation T which 
is dual to the Feynman diagram, while the edges of the boundary tetrahedra are labeled 
with the irreps Ai, Ag and A'^, Ag. 

The amplitude A(T) can be recognized as the Ponzano-Regge amplitude for the 
transition from an initial surface to a final one via triangulation T of the spacetime 
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manifold whose boundaries are the initial and the final surface, where the edges of the 
initial and the final triangulations are colored with a fixed set of irreps [1^ . 
The general spin net amplitude 

A,;(n) = ($(7',A')l(5.r|$(7,A)) , (77) 

can be computed from 

Af{n) = {<l>{^',A'){SXHl,^)), (78) 



by using the formula (|74|). It will be given by the sum of all (t>, v') FD with n ver- 
tices and V + v' external legs, multiplied by the normalization factor A^^^(A)A^^^(A'), 
where v and v' are the numbers of vertices in the spin nets 7 and 7', respectively. After 
summing over the representation indices, each FD will give the Ponzano-Regge am- 
plitude for the corresponding 3-complex T consisting of n tetrahedrons and with two 
disjoint boundaries consisting of v and v' tetrahedrons. One can then write 

Af{n)= J2 ^(^;7,A;7',A'). (79) 

r,\r\=n 

The colored boundary complexes Ty and 7J are dual to the spin nets 7 and 7' since 
they can be considered as two-complexes corresponding to triangulations of the initial 
and the final surface. 

Also, according to our time variable interpretation, the amplitude ( [79| ) will cor- 
respond to the transition from the initial to the final state in n units of the covariant 
time. 

Note that the spin network states |$(7, A)) have the same labeling as the spin net- 
work states |\l/(7, A)) obtained by the canonical quantization of 3d gravity, or equiv- 
alently, the SU{2) BF theory [[TT|]. However, these states are fundamentally different, 
because they belong to different Hilbert states, which is the consequence of two differ- 
ent quantization procedures which are used. 

In the canonical approach, the spacetime manifold M has a fixed topology R x S, 
where S is a 2d surface. In the connection representation, one can construct a Hilbert 
space of gauge invariant states L2{A/Q), whose orthonormal basis is given by the 
spin network states a(^)- These are evaluated by assigning the open holonomies 
U^{e, A) to the edges e of the spin net graph 7 which is embedded in S, where A is 
the irrep assigned to the edge e. One then takes the product of the holonomy matrices 
and contracts the indices with the intertwiners in order to obtain a scalar wave-function 

This construction can be discretized [|TU[ [TT]], where S is replaced by a triangulation 
T(E), and 7 is the dual graph, while the connection is replaced by a set of group 
elements g^, representing the holonomies of the edges e, so that U^p{e, A) — > D^p^g^). 
Then \E' becomes an element of the space L2{G^) where E is the number of the edges 



of 7 [|T4p. Hence in this formulation, one can consider only the spin nets with a fixed 
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number of the edges, so that one cannot calculate the transition amplitude between 
the spin nets having different numbers of edges. On the other hand, in the case of 
field theory spin networks, there is no such a constraint, because the Hilbert space has 
the Fock space structure, and one can have spin net states with different numbers of 
the edges. Hence the field theory Hilbert space contains states describing all colored 
2-complexes, or equivalently, the states for manifolds E of all topologies. 

Note that a single complex amplitude in (^9]) can be written as a scalar product of 
the canonical quantization spin net states |\E'(7, A)) as 

^(r,7,A;7',A') = (^(7',A')I^(7,A)), (80) 
when 7 and 7' have the same numbers of edges [[T5|, [T^. On the other hand, it follows 



from (79) that 



Af{n) = C(n,7,7')(^(7',A')l^(7,A)) 

= ($(7',A')|(5,r|$(7,A)) , (81) 

where C{n, 7, 7') is the number of 3-complexes T with n tetrahedrons whose bound- 
ary complexes are dual to 7 and 7', or equivalently, it is the number of (v, v') FD with n 
vertices. Hence due to the topological nature of 3d gravity, one can relate the canonical 
and the field theory spin net states. 



7 D=4 model 

The case of four spacetime dimensions is special for the BF theory, since this is the 
dimension where the BF theory starts to be different from the GR theory. GR iwD > A 
is not a topological theory, or in other words, it has local degrees of freedom. This is 
also reflected by the fact that the two-form Bab, where B = B"-^Jab and J^b are the 
5*0(4) or 50(3, 1) Lie algebra generators, cannot be always written as e"Ae^, where e" 
is the fierbein one-form. Recall that in 3d B"-^ = e"-^^ec, so that the BF action becomes 
identical to the Palatini form of the Einstein-Hilbert action. Therefore in order to 
obtain a simplical non-topological gravity theory from the BF theory, something has 
to be modified in a 4d analog of the constructions described in the previous sections. 

The field theory which generates the partition function of the topological gravity 
theory in 4d was worked out by Ooguri [^. The labeling pattern in the 4d action is 
a straightforward generalization of the D = 2,3 theories. We assign (f{gi, 5^4) to 
a tetrahedron whose faces are labeled with the group elements gi, and the interaction 
term is given by a 4- simplex built out of 5 tetrahedrons, so that 

Sa = \ I A<^'(1,2,3,4) (82) 

+ A /■ d'Og ^(1^ 2, 3, 4)<^(4, 5, 6, 7M7, 3, 8, 9)<^(9, 6, 2, 10)(^(10, 8, 5, 1). 
5! Jgw 
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Since Lp has to be G-invariant, we have 

'^{gi,g2,g3,g'i) = dg(f){gig,g2g,g3g,g4g) = PG(t>{gi,g2,g3,g4) ■ (83) 
Jg 

This gives the following momentum mode expansion 



^=J2 V^(Ai)ci(A2)rf(A3)rf(A4)0S^'<m4T'(^i'^?2, ^3,^74) , (84) 

A,A',a 

where 

, A(1234)A' _ , A(1234) ^A(1234)A' .o^x 

*^a(1234) ~ 2-^ *^a/3(1234) ^ f3{1234.) ' ^OJ; 

and 

/3 i=l 

The coefficients C are components of the intertwinwer map L4 

^4 : ^(Ai) ® V^(A2) ® ^(A3) ® V^(A4) ^ V^(O) , (87) 
in an orthonormal basis. 



In order to obtain ( |84| ) one needs the relation 

/ dg t/^^ {g) ® ^7^^ ((7) ® ?7^« ((?) ® ?7^^ ((7) = , (88) 

which in an orthonormal basis |a) becomes 

dgDtl,^ig)Dtl,M<,si9^^^^^^^^ > (89) 



G 



where A' labels the singlets in the tensor product of four irreps Ai, A4. 
By inserting ( ^4| ) into ^'4 we obtain for the kinetic part 



,A(1234)A' ^ 



2 |^a(1234) 
A,A'.Q 



where the 15j symbol is defined as 



J Ai ■ • • Aio 1 _ ^A(1234)Ai ^A(4567)A^ ^A(7389)A;j ^A(96210)A^ ^A(10851)A[; 
1 A' ■ ■ ■ A' f ~ ^a(1234) ^a(4567) ^a(7389) ^q:(96210) ^a(10851) ' 



(90) 



while for the interaction part we get 

q _ ^ Sr^ ,A(1234)A'i ,A(4567)Ai, ,A(7389)A^ , A(96210)A^ , A(10851)AJ; J Ai ■ ■ ■ Aio 
'-'^ ~ ^ *i^Q(1234) ^a{4567) ^a{73S°i) *^q:(96210) *i^a(10851) ] y\' ... \i 



(91) 



(92) 
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We have written the formulas ( Plj ) and for real D{g) matrices. In the complex 
case the repeated a indices will have a minus sign and there will be the sign factors 
coming from the complex conjugation rules for the D(g) matrix elements. Since this 
will not be important for our purposes, we will use the simpler real basis notation. 
The kinetic action (Q) implies for the propagator 



1=1 



Note that for a tetrahedron with labeled faces it is natural to consider ip which is in- 
variant under the permutations of the labels. This would then give the constraints on 
the Fourier modes (the analog of the hermiticity in 2d case), which would involve the 
6j symbols [H]. As a result, the propagator will not have the delta function for the 
intertwiner labels, but it will have the 6j symbol for Aj, A and A'. Since it is not clear 
whether the symmetric cp is necessary, and in order to make the presentation simpler, 
we will not impose the permutation symmetry. 

The Feynman diagrams will be given by the 5-valent graphs. These can be repre- 
sented as ribbon graphs, where each edge is replaced with four parallel lines, and each 
vertex is replaced with a pentagram whose comers are cut, and the incident lines for a 
corner are joined with the four lines of an edge. The value of a FD can be obtained by 
coloring the ten lines of each pentagram with ten irreps A, as well as every corner with 
the intertwiner label A'. One then assigns to each pentagram the corresponding 15 j 
symbol. For each propagator line one takes a delta function of the labels it connects. 
One then takes the product of all delta functions with the 15j symbols, and sums over 
the labels. 

Consider a vacuum FD T with \ V\ vertices, which is a dual one complex to a 4d 
simplical complex T consisting of \ V\ 4-simplices. For example, the hexagon graph, a 
six vertex 5-valent graph with all the vertices connected, is a dual one-complex for a 5- 
simplex considered as a 4d simplical complex. This 4d complex can be considered as a 
triangulation T of S"', which consists of five 4-simplices. Since each closed propagator 
line gives a d(A) factor, the Feynman rules associate to each loop (face) F of F an irrep 
Ap with the weight d{AF). For each vertex there is the 15j symbol of 10 loops (faces) 
which share it plus 5 intertwiner labels of the 5 edges coming out of the vertex. Hence 
one obtains 

z(r)= E n ''(M n I a''"' ::: 

Since a triangle / of T is dual to a face F of F, a tetrahedron t to the edge E, and 
a four-simplex s to a vertex V , we can rewrite the above expression as 

z(M) = E n ''(A/) n { If ::: a;";';'}- 
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Hence Z{M) is a sum over colorings of the complex T{M), where each triangle / of T 
is colored with an irrep Aj, while each tetrahedron t of T is colored with an intertwiner 
label from the tensor product of four irreps coloring the triangles of t. This sum is 
divergent, and by going to the quantum group formalism one obtains a finite number 
which is independent of the triangulation, and hence represents a topological invariant 
of the 4-manifold [ffl]. 



The form ( |95| ) of the partition function is the state sum which mathematicians have 
explored over the years, and it was in this approach that the modification which lead 
to quantum GR was discovered [|]]. Namely, the constraint B'^^ = e° A e^, which 
transforms the BF action into Pallatini action, can be written in a simple form as 

B'''aB,, = , (96) 



where the indices are contracted with the group metric. The constraint ( |96[ ) is translated 
in the state-sum formalism into a constraint on the representations with which one 
colors the triangles of a triangulation as 



r'jab = , (97) 

where J"^ are the generators of the 50(4) Lie algebra so(4). Since so(4) = su{2) © 



su{2), the constraint ( |97| ) becomes 



J^~K^ = , (98) 

where and K"^ are the Casimirs of the two SU{2). 

In terms of representations this constraint implies that one labels the triangles / 
of a 4-simplex s with the simple irreps = (j, j), and instead of summing over all 
irreps A = (j, k) in the state sum, one sums only over the simple irreps. Then the 15j 
symbol becomes a sort of lOj symbol, the Barret-Crane vertex 

^)^ AT AT \ — \^ ^Af{1234)Mi ^Af(4567)M2 ^Af{7389)M3 ^Ar{96210)A/4 ^Af(10851)M5 
|/(^iVi • ■ ■ iVioJ — 0^(^234) *-"a(4567) *-'a(7389) *-'a(96210) ^ a{W85l) ' 
M,a 

(99) 

where M labels the simple irrep singlets appearing in the tensor product of four simple 
irreps. Hence 

Nf fer seT 

where we have written Z as a function of T, since now the state sum will depend on 
the local parameters of the triangulation T. 

An important property of the BC vertex is that it can be expressed in terms of an 
integral over SU{2f Q 

v(iVi---iVio) = / d'h rr Tvj^^ihkh;') , (101) 

JSU{2)^ l<fc<K5 
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where (ji2, j45) = (A^i, A^io) label the edges of a 4-simplex s, and hi label the 
vertices. This formula can be easily proved by noticing that the right-hand expression 
is [10j(5f/(2))]2 and since 50(4) ~ SU{2) x SU{2), then lOj{SO{A))\j=simpie = 
[10j(5f/(2))]2. 

The explanation of this formula comes from the fact that 

SU{2) ~ = S0{A)/S0{3) , 

so that the h variables are really coordinates x on the homogeneous space X = G/H, 
which is a consequence of the fact that = (j, j) is a class-one representation of 
G = 50(4) with respect to the subgroup H = 50(3) 0. The quantity Trj{hkh-^) 
is the zonal spherical function, and can be thought of as a Green's function on the X 
space, so that the BC vertex can be represented as a Feynman graph of a field theory 
which lives on X. One then obtains 

Gn{xuX2) = Trjihih^^) (102) 

= ^m±})M , (103) 

sm Oi2 

where 9i2 is the angle between the unit 4-vectors Xi, or equivalently, the geodesic 
distance between the points of X. The vectors Xi can be interpreted as the normals to 
the tetrahedrons U which share the face /12 labeled with the irrep j. 

The representation ( |101D combined with the formula ( |103| ) gives the asymptotics 
of V in the limit of large j 



V 



5^A(A)cos(5^A;ey + 7r/4] , (104) 

A \f£s J 



where Aj = 2jf + 1 is the area of the face /, and one has to sum over certain in- 
equivalent colorings A of the four-simplex s, for details see [[TP|]. When ( |104| ) is in- 



serted in Zbc^ it is easy to see that the terms of the form Qe'^^ will appear, where 
Sr = J2fer ^f^f i^ the Regge form of the EH action. This is an indication that the 
semi-classical limit of the BC model is a theory that is related to classical GR. How- 
ever, the semi-classical limit of the BC model has to be investigated in more detail in 
order to be able to make more definite statements. 

Another important consequence of the formula ( |10ip is that it can be easily gener- 
alized to the Lorentzian case [0], since then G = 50(3, 1) ~ 5L(2, G), H = 50(3) 
and therefore X is a hyperboloid. The unitary irreps of 5L(2, G) are labeled with two 
numbers, i.e. A = {n,p) where n E Z and p G R, and the simple representations are 
(0, p) and (n, 0). It is natural to assign (0, p) irreps to space-like triangles and (n, 0) to 
time-like triangles. One can restrict the model to the continuous irreps only so that 

sinipdu) 

Gp{xi,X2) = — ■ , , , (105) 
p smn 0,12 
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where is the geodesic distance between xi and X2. By using ( JlOSp one can construct 
the Lorentzian analog of the BC vertex via 

V{N,■■■N^o)= f (fx J] Gp^X^u.xi) , (106) 

■^^^ l<k<l<5 

where (A^i, iVio) = (j9i2, •••,^45)- 

As discussed in [ [T8| ] for the Euclidean case and in for the Lorentzian case, 
the constructions ( \101\ ) and ( |106[ ) can be generalized to give expressions for higher 



nj symbols. Given a graph 7 we can label the vertices of 7 with the coordinates Xk 
from the homogeneous space X = G/H. Then associate to each edge (fc, Z) of 7 a 
simple irrep labeled with jki- In this way one obtains a labeled graph which is called 
a relativistic spin net. We can than associate the "propagator" Gj^., (xjt, a;;) to each 
edge, and if we now consider 7 as a Feynman diagram of a theory with delta function 
vertices, then its value will be given by 

V{N,---Nr,)= f (Tx J] G,^X^k.xi) , (107) 

l<k<l<m 

where < k < I < m} = {Ni, ■ ■ ■ , N^}. 

In the Euclidean case, the expression ( |107D is finite, since it is equal to a finite sum 
of products of a finite number of C symbols. This is a consequence of the formula 
( |102| ). In the Lorentzian case, the finiteness is not obvious. The lOj symbol was 
conjectured to be finite in [^, and recently a proof of finiteness was given in p^]. 

As far as the finiteness of the state sum Zbc{'^) is concerned, it is divergent, but 
as argued in one expects to become finite in the quantum group formalism. It is 
interesting that there is a modification of the BC field theory model in the Euclidean 
case which seems to give finite Z for all T . 

8 BC field theory 

Now we can discuss the field theory formulation of the BC model [|12|]. The constraint 
of using only the simple representations to label the triangles of a triangulation is 
translated in the field theory formalism into a requirement of invariance under the 
5*0(3) subgroup: ^{gi) = ^{gihi) , hi e H. Hence 

^{91,92,93,94) = d^h(f){gihi,g2h2,g3h3,g4h4) = Ph4> , (108) 

where (p is an unconstrained field. By combining the iJ-invariance ( |108| ) with the usual 
G-invariance (BSb one obtains 



(p = PgPh<P = dg d^h 4>{gighi) 
Jg Jm 

= y: vkn.) ■ ■ ■ 4iV4)<s' n . <io9) 



N,a,l3 1=1 
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where 



B 



iV(1234) 



EC, 



Ar(1234)M 
/3(1234) 



(110) 



M 



M denotes the simple irreps whose singlets appear in the tensor product of four simple 
irreps, and A is the dimension of that subspace. By inserting ( |109| ) into 5*4 one obtains 



Sk 



N,a 



,Af(1234) 
a(1234) 



(111) 



while for the interaction part we get 



A 



^ ^^(1234) , iV{4567) , 7V(7389) , 7V(96210) , 7V(10851) ^ ^ x ... 

^ V^a(1234) '/^«(4567) V^a(7389) V>a(96210) V>a(10851) "^l^^l ' ' ' ^^loj , (H^) 



where V is the normalized BC vertex 



V 



V(iVi---iVio) 



v/A(1234)---A(10851) 



(113) 



Hence the FD of the field theory defined by ( |112| ) will be given by the closed and 
open 5-valent graphs, as in the topological case, but the propagator and the vertex 
factors will be different. A vacuum FD T will reproduce the state sum ( |100D for a 
triangulation T whose dual one-complex is T. However, now one expects that the 
regularized sum will depend on the triangulation T, since Z{T) does not have the 
properties of a topological invariant. 

Note that there exists a modification of the BC field theory action which gives finite 

1. The 



FD in the Euclidean case without using the quantum group as a regulator [ pl| , | 
modified action has the same kinetic term as the topological action (G invariance only), 
while the interaction term is the same as in the BC case (G and H invariance). When 
expanded in the Fourier modes, one obtains 



NM,a 



{dN^...di^^Y 



, Ar(1234)Af 
o(1234) 



(114) 



for the kinetic term, while 



5! 



E 

N,M,a 



,Af(1234)Afi 
«(1234) 



,Af(4567)Af2 
a(4567) 



,Af(7389)Af3 
a(7389) 



,Af(96210)Af4 
a(96210) 



7V(10851)A/5 
a(10851) 



10) 



is the interaction term. The appearance of the factor {d]\f-^...d 



N4 



(115) 

)^ in the kinetic term 
gives a damping factor (dN^...diy^)^'^ in the propagator which insures the finiteness. 

As far as the FD with external legs are concerned, they will correspond to trian- 
gulations of a 4-manifold with boundaries, where each boundary contains more then 
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one 4-simplex. The boundary complexes are then interpreted as 3d complexes, corre- 
sponding to triangulations of boundary 3 -manifolds Ej. These FD will be relevant for 
the transition amplitude from a state on a 3-manifold Si to a state on a 3-manifold E2, 
where dM = Si U S2. However, we need to construct first the spin net state associated 
to a triangulation of a 3-manifold S representing a boundary of M. This will be given 
by a straightforward generalization of the results from the lower dimensions. 

Consider a 4-valent connected closed graph 7 which corresponds to a triangulation 
of S, i.e. 7 is a dual one complex to T(S), where now T(S) is considered as a 
3d simplical complex. For example, a 4-simplex is dual to a pentagram. A more 
complicated example is a 3d simplical complex consisting of 8 tetrahedrons. The dual 
graph will have eight 4-valent vertices and 16 edges. 

Given a triangulation T(S) we can construct a functional $ by integrating a prod- 
uct of ^i{gj) over G", where n is the number of triangles in T(S) and ^i{gj) is asso- 
ciated to the tetrahedron in T(S) whose triangles are colored with the group elements 
9j, i.e. 

^ = / d^^g n v{9Mt),gf2{t),9fB{t),9Mt)) ■ (116) 
By Fourier expanding it, we get 

^ = E • • • C7:;r>7(A^iA^2 ■ ■ ■ at^j , d n) 

where the vertex V-y(A'^i ■ ■ ■ A^„^) is given by ( |1U7D and 7 is the four-valent graph asso- 
ciated to T(S), i.e. a dual one-complex to T(S). Hence we take the operator 

$(7,iV) = 0^(1234^^ 

~ jVnfA^) ^-'■l'^a(ei(f)...e4M) '^^^^^ 

to represent the spin net given by the graph 7 whose edges are labeled with the simple 
irreps A'^i, Nn^- In the case of the Perez-Rovelli model, one obtains 

$(^, AT) = y rr ^A^(eiM...e4(.))M„ ^ ^ 

The amplitude 

Ai2H = ($(72,iV2)(5.)"$(7i,iVi)) , (120) 

will be interpreted as the amplitude for a transition from the state |$(7i, A'^i)) to the 
state [$(72, ^^"2)) in n units of the covariant time T. This will be a finite sum of FD, 
consisting of {v\,V2) FD with n vertices and v\ + V2 external legs, where Vi is the 
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number of the vertices in the spin net 7j. Each FD would give a 4d analog of the 
single-complex Ponzano-Regge amplitude, if we choose a normalization Nq(N) = 

Note that the amplitude ( |120| ) will be finite in the Euclidean Perez-Rovelli model. 
Therefore this would be en example of a finite 4d quantum gravity theory with local 
degrees of freedom, provided that one uses the covariant time interpretation. In the 
Lorentzian case it is not known whether the Perez-Rovelli model is finite, but it is 
plausible that it may be finite. Alternatively, one can try the quantum Lorentz group 
regularization. 

As in the 3d case, the field theory spin nets $ can be considered as generaliza- 
tions of the canonical quantization spin nets \E' [ p4| ] to situations when the topology 
of the spacetime manifold is not fixed. Also, due to the non-topological nature of 4d 
gravity, the labels are not the same anymore. Covariant spin nets have simple 5*0(4) or 
SL(2, C) irreps as labels, while the canonical spin nets have the SU (2) labels. In order 
to make a closer link between these, one would have to discretize the canonical spin 
nets, as in the 3d case, but in the 4d case one would have to solve the diffeomorphism 
and the Hamiltonian constraints. Furthermore, one would have to find a time variable 
in order to compute the transition amplitudes. A great advantage of the covariant ap- 
proach, which is related to the fact that its a path-integral quantization, is that one does 
not have to do all these non-trivial steps, and one can calculate directly the transition 
amplitudes. 



9 Conclusions 

We have demonstrated that there is a significant benefit if one completes the Feyn- 
man diagram picture of the transition amplitudes for the state- sum models of quantum 
gravity by introducing the Fourier mode operators and the corresponding Fock space, 
so that the transition amplitudes can be understood as matrix elements of an evolu- 
tion operator, in exact analogy with the particle field theory case. The spin network 
operators are constructed from the products of the field Fourier modes, such that each 
mode corresponds to a vertex of the spin net graph, and then one averages over the spin 
components in a group-theory invariant way. The corresponding state can be thought 
of as a state of a universe whose space is discrete. This 3d space is given by the 
dual 3-complex to the spin net graph 7, and the quantum states are determined by the 
colorings of the triangles with the irreps A. 

We have considered only the transition amplitudes for the single-universe states. 
Since one can act on the Fock space vacuum with more than one spin net operator, 
one can obtain the multi-universe states, and the corresponding transition amplitudes. 
Hence our quantum field theory of spin networks is also a quantum field theory of 
discrete universes. This QFT is then a discrete realization of the third quantization of 
gravity. 

Note that the expressions for the amplitudes we have written are mostly formal. 
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and they need a regularization for the infinite sums over the irreps of G. The standard 
approach is to use the quantum group formalism, which is a group invariant way of 
putting a cut off on A. Alternatively, the Euclidean BC model can be regularized via 
Perez-Rovelli reformulation. Hence in this model the amplitudes for finite time inter- 
vals will be finite. The same could apply for the Lorentzian version of the Perez-Roveli 
model [p2|], provided one proves the finiteness of an arbitrary Feynman diagram. If this 
turns out to be true, one would have a well defined quantum theory of 4d Lorentzian 
gravity. Then the main task would be to examine the semi-classical states of this the- 
ory. Given the discrete nature of the theory, one expects that such states would be 
described by spin nets of large number of vertices, since large triangulations represent 
good approximations of smooth manifolds. 

The key ingredient in making the quantum field theory of spin nets well defined 
is the introduction of the discrete time variable T and the corresponding evolution 
operator U (T). Since U (T) = (S^)'^ one can define the transition amplitudes for finite 
time intervals, which are given by the sum of FD with T vertices, or equivalently, by 
the sum of partition functions for triangulations with T simplices. In this way one 
avoids the infinite sum of the single-complex amplitudes, or Feynman diagrams, for 
all possible complexes interpolating between the boundary complexes, which is often 
invoked in the literature as the physical amplitude. We have demonstrated that this is 
not necessary, since this prescription is related to the evolution operator U = exp^iSv), 
and this operator is not natural in the covariant quantization formulation of theories 
without a fixed spacetime background. 

The construction of the covariant time evolution operator mirrors the physical in- 
tuition that summing over all possible triangulations of the spacetime manifold corre- 
sponds to including metrics which give arbitrary large geodesic distance between the 
initial and the final boundary. Therefore one should divide this sum into parts corre- 
sponding to fixed time-like geodesic lengths. These partial sums still contain infinitely 
many triangulations. Hence a better time variable, i.e. a variable which gives finite 
partial sums, is the discretized analog of the spacetime volume, which is the num- 
ber of the Z^-simplices in a simplical complex representing the spacetime. Also, this 
interpretation is natural for the compact universes. 

In order to gain a better understanding of the covariant time, one would have to see 
what is its relation to canonical time variables[], as well as a relation to matter clock 
variables. This analysis will be also important for understanding the meaning of the 
non-unitarity of the covariant time evolution operator. 

If one takes the covariant evolution operator as the fundamental object, this would 
mean abandoning the law of unitary evolution of states in quantum mechanics (QM). 
However, this may not be a bad thing in this context, since we are dealing with the 
wavefunction of the Universe, which as a concept is not well defined in the standard 

^These are time variables in the canonical formulation of GR [Q. In this context a typical canonical 
time variable would be the volume of a spatial section, which would correspond to the number of vertices 
in the spin net. 
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QM due to the fact that there is no outside observer. Also, it could provide a mecha- 
nism for the the matter wavefunction collapse, which would be a concrete realization 
of gravity driven wavefunction collapse Also note that the QFT formulation of 
spin nets is also a good arena for exploring other interpretation schemes of QM for 
quantum gravity [P5|]. 

In order to explore all these ideas, as well as in order to have a complete theory, 
one would need a simplical field theory formulation of matter coupled to gravity. 
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